NPS ARCHIVE 
1997.06 
BUCH, E. 

NPS-OC-97-004 

NAVAL POSTGRADUATE SCHOOL 
MONTEREY, CALIFORNIA 




THESIS 



WIND-FORCED MODELING STUDIES OF 
CURRENTS, MEANDERS, EDDIES, AND 
FILAMENTS OF THE CANARY CURRENT 
SYSTEM 
by 

Eric J. Buch 
June 1997 

Thesis Advisor: Mary L. Batteen 

Approved for public release; distribution is unlimited. 



Thesis 

B82944 



Prepared for; 

’■’"'tional Science Foundation 
(1 Wilson Boulevard 
ington, VA 22230 



'yUi.EVK'NO/ .iBfWtY 
pAi. fOSTGRADUATESCHOO 
lUAlTEREN Of W943-5101 



NAVAL POSTGRADUATE SCHOOL 



ooDLEVKNOXUBRARy MONTEREY, CALIFORNIA 

NAVAL POSTGRADUATE SCHOOL 
MONTEREY, CA 939434101 



Rear Admiral Marsha J. Evans 
Superintendent 



This thesis was prepared in conjunction with research sponsored in part by the 
National Science Foundation, 4201 Wilson Boulevard, Arlington, VA 22230. 

Reproduction of all or part of this report is authorized. 



Released by: 



REPORT DOCUMENTATION PAGE 



Form Approved 0MB No. 0704-01 8S 



Public reporting burden for this collection of mformatioQ is estimated to average 1 hour per response, including the tim e for reviewing instruction, searching existing data 
sources, gathering and m ai n ta in i n g the data needed, and completing and reviewing the collection of information. Send comments regarding thLo burden estimate or any 
other asfwv-l of this collection of information iochkling sugge^tiunh redurin§ iluj hkudtz.. to Wulingtjii Kcndx)uJitt«r<r Services, Dva^iotbIb for biAmnslioc 
Operations and Reports, 1215 Jefferson Davis Highway, Suite 1204, Arlington, VA 22202-4302, and to the Office of Management and Budget, Paperwork Reduction 
Project (0704-0188) Washington DC 20503. 



1 . AGENCY USE ONLY (Leave blank) 


2. REPORT DATE 3. REPORT TYPE AND DATES COVERED 

June 1997 Master’s Thesis 


4. TITLE AND SUBTITLE WIND-FORCED MODELING STUDIES OF 
CURRENTS, MEANDERS, EDDIES, AND FILAMENTS OF THE 
CANARY CURRENT SYSTEM 


5. FUNDING NUMBERS 


6. AUTHOR(S) Eric J. Buch 


7. PERFORMING ORGANIZATION NAME(S) AND ADDRESS(ES) 
Naval Postgraduate School 
Monterey, CA 93943-5000 


8. PERFORMING 

ORGANIZATION 
REPORT NUMBER 
NPS-OC-96-004 


9. SPONSORING/MONITORING AGENCY NAME(S) AND ADDRESS(ES) 
National Science Foundation 
4201 Wilson Boulevard Arlington, VA. 22230 


10. SPONSORING/MONITORING 
AGENCY REPORT NUMBER 


11. SUPPLEMENTARY NOTES The views expressed in this thesis are those o1 
official policy or position of the Department of Defense or the U.S. Govern 


F the author and do not reflect the 
ment. 


12a. DISTRIBUTION/AVAILABILITY STATEMENT 

Approved for public release; distribution is unlimited. 


12b. DISTRIBUTION CODE 


! 13. ABSTRACT (maximum 200 words) 

1 A high-resolution, multi-level, primitive equation ocean model is used to examine the response 

■of an eastern boundary oceanic regime to both wind forcing and irregular coastline geometry. The 
focus of this study is the coastal region from 30° N to 42.5° N, a portion of the Canary Current System 
(CCS). To study the generation, evolution, and sustainment of the currents, meanders, eddies and 
filaments of the CCS, the model is forced from rest using seasonal climatological winds. To investigate 
i the role of irregular coastline geometry, the first experiment uses climatological wind forcing along an 
idealized “straightened” coastline, while the second experiment uses the same wind forcing along an 
irregular coastline. In both cases a surface current, undercurrent, meanders, eddies, and filaments are 
; generated. The results obtained while using the irregular, rather than the idealized coastline, however, 
show preferred eddy generation locations as well as enhanced growth of meanders, eddies, and 
filaments. The features produced by the model are consistent with available observations of the CCS. 
The model results support the hypothesis that both wind forcing and irregular coastline geometry are 
1 important mechanisms in the generation of many of the observed features of the CCS 




SUBJECT TERMS Primitive equation model, Canary Current System, currents, 

meanders, eddies, filaments 


15. NUMBER OF 

PAGES 123 

16. PRICE CODE ! 


17. SECURITY CLASSIFICA- 
TION OF REPORT 
Unclassified 


18. SECURITY CLASSIFI- 

CATION OF THIS PAGE 
Unclassified 


19. SECURITY CLASSIFICA- 
TION OF ABSTRACT 
Unclassified 


20. LIMITATION OF i 
ABSTRACT 
UL 



NSN 7540-01-280-5500 



Standard Form 298 (Rev. 2-89) 
Prescribed by ANSI Std. 239-18 298-102 



11 



Approved for public release; distribution is unlimited. 



WIND-FORCED MODELING STUDIES OF CURRENTS, MEANDERS, 
EDDIES, AND FILAMENTS OF THE CANARY CURRENT SYSTEM 



Eric J. Buch 

Lieutenant, Uni'/ed States Navy 
B.S., United States Naval Academy - 1988 

Submitted in partial fulfillment 
of the requirements for the degree of 

MASTER OF SCIENCE IN METEOROLOGY AND PHYSICAL 

OCEANOGRAPHY 

from the 

NAVAL POSTGRADUATE SCHOOL 
June 1997 



DUDLEY KNOX LIBRARY 
NAVAL POSTGRADUATE SCHOOL 
MONTEREY, CA 93943-5101 



ABSTRACT 



A high-resolution, multi-level, primitive equation ocean model is used to examine 
the response of an eastern boundary oceanic regime to both wind forcing and irregular 
coastline geometry. The focus of this study is the coastal region from 30° N to 42.5° N, a 
portion of the Canary Current System (CCS). To study the generation, evolution, and 
sustainment of the currents, meanders, eddies and filaments of the CCS, the model is 
forced from rest using seasonal climatological winds. To investigate the role of irregular 
coastline geometry, the first experiment uses climatological wind forcing along an 
idealized “straightened” coastline, while the second experiment uses the same wind forcing 
along an irregular coastline. In both cases a surface current, undercurrent, meanders, 
eddies, and filaments are generated. The results obtained while using the irregular, rather 
than the idealized coastline, however, show preferred eddy generation locations as well as 
enhanced growth of meanders, eddies, and filaments. The features produced by the model 
are consistent with available observations of the CCS. The model results support the 
hypothesis that both wind forcing and irregular coastline geometry are important 
mechanisms in the generation of many of the observed features of the CCS. 
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I. mXRODUCTION 



The eastern boundary current (EEC) system off the west coast of the Iberian 
Peninsula and Northwest Africa is part of the Canary coastal upwelling system, which 
extends from -10° N to 45° N. There is general agreement (e.g., Wooster et ah, 1976^ 
that the Canary Current is generated as the closing eastern boundary for the North 
Atlantic Gyre, enhanced with northerly winds during and after the Spring Transition and 
reduced as the northerly winds relax or become westerly in the winter season or milder 
Fall Transition. The climatological mean Canary Current system (CCS) is a classical EEC, 
which consists of several large-scale currents. The predominant flow is the Canary 
Current (CC), which like other EECs is a broad (-1000 km), relatively slow (-10-30 
cm/s), equatorward, yearlong surface flow extending to depths of - 500 m (Wooster et 
ah, 1976). The extension of the CC which lies off the western coast of the Iberian 
Peninsula (BP) is also occasionally referred to as the Portugal Current (e.g., Tomczak and 
Godfrey, 1994). Embedded in the CC are narrower, somewhat swifter flows (occasionally 
reaching speeds of 100 cm/s (DMAH/TC, 1988)). Eelow the CC, lies the poleward 
undercurrent, which is narrower (-10-40 km) and weaker (-2-10 cm/s) than the surface 
current, and found near the coast. The depth and strength of the undercurrent varies both 
seasonally and latitudinally but is usually strongest near 300 m. In the winter the 
undercurrent shoals to the north, occasionally reaching the surface to form a third flow 
component. This current, which has recently been referred to as the Iberian Current (IC) 
(Haynes and Earton, 1990) is a seasonal surface poleward current, found normally north 
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of Cabo da Roca (see Figure 1 for geographical reference) but occasionally extending 
south nearly to Cabo de Sao Vincent during the winter months. The IC is also narrow 
(25-40 km) and relatively weak (-20-30 cm/s) and is found trapped near the coast against 
the shelf break (Fiuza, 1980; Frouin et al., 1990; Haynes and Barton, 1990). A second 
deeper undercurrent also flows poleward off the western coast of the Iberian Peninsula. 
This flow is attributed to the Mediterranean Outflow (MO) which has at least two main 
cores, a shallow core at depths of 600-900 m and a deeper core at 1100-1200 m (Ambar, 
1980). A third shallower poleward core of Mediterranean water has also been shown to 
exist around 400 m (Fiuza, 1979; Ambar, 1980). 

The CCS is influenced predominantly by equatorward, upwelling favorable winds 
produced by the Azores High. The Azores High is a semi-permanent subtropical high 
pressure system that is similar in nature and behavior to its counterpart off the west coast 
of North America, the North Pacific Subtropical High, as described in Nelson (1977). The 
center of the Azores High migrates meridionally with the seasons, reaching its 
southernmost extent, -27® N in March and ridging north to -33® N by August. Because 
of this migration, maximum values of wind stress vary temporally at given locations. 
During the summer the mean east-west pressure contrast between Portugal and the center 
of the Azores High is -8 mb, while in the winter it weakens to -1 mb. This results in 
considerably stronger northerly to northwesterly winds in the summer followed by weaker 
northwesterly winds to winds with even a slight southerly component off of northern 
Portugal and northwest Spain in the winter. The shift of maximum wind stress also causes 
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the upwelling favorable winds to shift from about 27° N near the Canary Islands in 
January, to about 43° N off the IP by July (Fiuza, 1982). 

Recent observations have shown that superimposed on the broad, climatological 
mean flow in this and other EBCs are highly energetic, mesoscale features such as jet-like 
surface currents, meanders, eddies, and filaments. Satellite images of sea surface 
temperatures have shown several filaments extending off the coast of Portugal and 
northwest Spain (Fiuza et al, 1989) as well as off Cape Ghir in northwest Africa (Van 
Camp et al, 1991; Hagen et al, 1996). Observations have also shown pairs of 
anticyclonic and cyclonic mesoscale eddies on the order of 100 km (Fiuza, 1984) off the 
coast of Iberia. These features have been observed during periods of predominantly 
equatorward winds, which are upwelling favorable and many seem to be located near 
prominent coastal features such as capes. These observations provide evidence that wind 
forcing and coastal irregularities are possibly important mechanisms for the formation and 
sustainment of many of the mesoscale features found in EBC regions. 

There have been numerous wind forcing models over the past few decades of 
EBCs, particularly for the California Current System. Early work included that of 
Pedlosky (1974) and Philander and Yoon (1982), who used steady wind stress and 
transient wind forcing, respectively. Carton (1984) and Carton and Philander (1984) 
investigated the response of reduced gravity models to realistic coastal winds. Another 
series of experiments were conducted by McCreary et al (1987) utilizing a linear model 
with both transient and steady wind forcing in the California Current System. In all these 



3 



previous models, relatively weak currents (5-10 cm/s) were generated and no eddies or 
filaments developed. 

Recent modeling studies have focused on the driving mechanisms for the complex 
features such as filaments observed in the EEC regions. While Ikeda et al. (1984) and 
Haidvogel etal. (1991) studied baroclinic and barotropic instability, coastal irregularities, 
and bottom topography as possible mechanisms, Batteen et al. (1989), McCreary et al. 
(1991), Pares-Sierra et al. (1993), and Batteen (1997) studied wind forcing as a possible 
generative mechanism. Recently, Batteen (1997) used a multilevel primitive equation (PE) 
model run with seasonal winds and coastal irregularities to realistically generate many of 
the observed features of the California Current System. 

In contrast, modeling studies of the CCS have been scarce. McClain et al. (1986) 
performed the first limited modeling study in the region, using ship-derived winds to 
produce a large negative wind stress curl off the northwest coast of the IP. This resulted 
in opposing equatorward and poleward surface currents. Mittlestaedt (1991) looked at 
the near-surface circulation of the PE models of Bryan (1969) and Semtner (1974). 
Because these coarse resolution models could only resolve structures larger than 200 km, 
they missed the mesoscale features. Recently Batteen et al. (1992) used an eddy resolving 
PE model with both uniform and variable wind stress, the latter computed from synoptic 
surface pressure analyses, to produce realistic currents as well as mesoscale eddies. This 
study, however, covered only a limited region of the CCS off the northwest coast of the 
Iberian Peninsula and also used a straight coast for its model domain. 
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The objective of this study is to investigate the dependence of the generation of 
currents, eddies, jets and filaments on seasonal climatological wind forcing and on 
irregular coastal geometry for an extended region of the CCS. The high-resolution, multi- 
level, PE model of Batteen (1997) with seasonal wind forcing will be used. As in the 
California Current systems study, to allow larger scale eddies and elongated filaments to 
be generated and so that the model can reach equilibrium, the model will be allowed to run 
for ~3 years, significantly longer than in previous studies which ran for less than a year. 

The study is organized as follows: The numerical model, types of forcing, initial 
and specific experimental conditions, and the energy analysis technique are described in 
section II. The results of the different model simulations and comparisons \vith 
observations are discussed in section III, while a summary and discussion are presented in 
section IV. 
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n. MODEL DESCRIPTION 



A. MODEL EQUATIONS 

The PE numerical model in this study was originally a coarse resolution model 
used in closed basin studies by Haney (1974). It has recently been adapted for eddy- 
resolving, limited EEC regions with open borders on the northern, western, and southern 
boundaries by Batteen (1997). The limited area EEC model is multi-level, non-adiabatic, 
and uses the beta-plane approximation. It has both baroclinic and barotropic velocity 
components and uses the hydrostatic, Boussinesq, and rigid lid approximations. The 
governing equations are as follows: 



du - \ ^ . d^u 


( 1 ) 


dv . d^v 


( 2 ) 


o 

II 


( 3 ) 


II 


( 4 ) 


P = Po[^-a{T-T^) + P{S-S^)] 


( 5 ) 


dt " " 


( 6 ) 


d, ■ ^ 4" 


( 7 ) 
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In the above equations, l is time, {x,y,z) is a right-handed Cartesian coordinate 
system with x pointing toward shore, y alongshore, and z upward. The corresponding 
velocity components are (w,v,w), T is temperature, S is salinity, p is density, and p is 
pressure. Table 1 provides a list of other symbols found in the model equations, as well as 
values of constants used throughout the study. 

A space-staggered B-scheme (Arakawa and Lamb, 1977) is used for the horizontal 
finite differencing. This scheme has been shown by Batteen and Han (1981) to be 
appropriate when the grid spacing is approximately on the same order as, or less than, the 
Rossby radius of deformation. The horizontal grid spacing is 9 km in the east-west 
direction and 11 km in the north-south direction, while the internal Rossby radius of 
deformation is ~30 km. The model uses ten layers in the vertical, with constant z-levels, 
at depths of 13, 46, 98, 182, 316, 529, 870, 1416, 2283, and 3656 m. This spacing 
scheme concentrates more layers above the thermocline in the dynamically active portion 
of the ocean, consistent with Haney (1974). 

The model domain (Figure 1) is rectangular, encompassing the west coast of the 
Iberian Peninsula and northwest Africa, from 30° N to 42.5° N (1408 km alongshore), and 
from 5° W to 17.5° W (1152 km cross-shore). The coastal boundaries of the model 
domain are closed, and have both the tangential and normal components of velocity set to 
zero. To focus on the roles played by climatological wind forcing, bottom topography has 
been omitted and the eastern boundary is modeled as a vertical wall. The constant depth 
used in the model is 4500 m. A modified version of the radiation boundary conditions of 
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Camerlengo and O'Brien (1980) is used for the open ocean domain boundaries to the 
north, south, and west with slight spatial smoothing applied in the vicinity of the open 
boundaries (Batteen, 1997). 

Biharmonic lateral heat and momentum diffusion is used in the model with the 
same choice of coefficients (i.e., 2.0 xio’’ cmVs) as in Batteen (1997). Holland (1978) 
showed that the highly selective biharmonic diffusion acts predominantly on 
submesoscales, while Holland and Batteen (1986) found that baroclinic mesoscale 
processes can be damped by Laplacian lateral heat diffusion. As a result, the use of 
biharmonic lateral diffusion should allow mesoscale eddy generation via barotropic 
(horizontal shear) and/or baroclinic (vertical shear) instability mechanisms. As in Batteen 
(1997), weak (0.5 cm^/s) vertical eddy viscosities and conductivities are used. Bottom 
stress is parameterized by a simplified quadratic drag law (Weatherly, 1972), as in Batteen 
(1997). 

An annual climatological temperature profile based on Levitus (1982) centered at 
~37.5° N (corresponding to the center of the model domain) was used to represent the 
initial temperature conditions for the model domain. An exponential profile was used to 
approximate the initial mean stratification in both experiments. The exact form was: 

T(z) = Tb + 
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The approximation uses a length scale of 450 m, and assumes Tb = 2° C to be the 
temperature at depth, and AT = 15.5 °C to be the change of temperature between the 
bottom and top of the profile. 

To isolate the effects of wind forcing in this experiment, the net heat flux at the sea 
surface is zero. Therefore, any heat flux that is discernible during the model runs will have 
been generated by an increase or decrease of sea surface temperature resulting from wind- 
forcing effects (Batteen et al., 1989). To accomplish this, an initial air temperature is 
chosen that forces the net flux of longwave radiation, sensible heat, and latent heat to 
balance the heating due to solar radiation. This air temperature is then used in the model 
for all wind forcing experiments. Any subsequent surface heat flux forcing is therefore a 
secondary effect of the changes to the sea surface temperature due to the wind forcing 

B, METHOD OF SOLUTION 



Equations (1) through (7) comprise a closed system of seven scalar equations and 
seven unknowns, u, v, w, p, p, T, and S. The variables, u, v, T, and S are prognostic 
variables whose time rates of change are predicted from (1), (3), (6) and (7), respectively. 
Although the diagnostic variables w, p, and p can be determined from (3), (4), and (5), 
respectively, additional constraints are imposed on p and w by the choice of the rigid lid 
boundary conditions. The vertically integrated pressure can no longer be obtained by 
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integrating the hydrostatic equation (4) for the free surface, and the vertically-integrated 
horizontal velocity is subsequently constrained to be non-divergent, i.e.. 



{du 



de = Q , 



( 8 ) 



which is obtained by integrating (3) and applying the vertical boundary conditions where 
£■ is a dummy variable representing the vertical coordinate. 

For any quantity q, let its vertical average be denoted by q and its departure 
(vertical shear) by q' . From (8) the vertical mean flow can then be described by a 



streamflmction y/, such that; 



u = - 



1 



_ _ \ dy/ 



( 9 ) 

( 10 ) 



The streamflinction is predicted from the vorticity equation, which is derived by 



applying the curl operator to the vertical average of (1) and (2), and then using (9) and 
(10), the vorticity equation becomes 
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where G and F represent the collected contributions of the nonlinear and viscous terms 
from equations (1) and (2). 

The vorticity equation (11) is solved by obtaining an updated value of ^ by 

application of the leapfrog (or every 1 1 time steps, the Euler-backward) time-differencing 
scheme. The associated value of y/ can then be obtained from: 






1 ( 1 r dH 

\ dx^ 



H\dy^ 



-1 



dx dx. dy dy 



( 12 ) 



which is an elliptic equation. A solution to (12) is fully prescribed by specifying the values 
of y/ on the open and closed boundaries of the model domain. Currently, to solve (12), 

the model uses an elliptic solver when there are no variations in coastline geometry and/or 
topography, and successive over-relaxation techniques when there are variations in 
coastline geometry and/or topography. 

The vertical shear current v') is predicted from (1) and (2) after subtracting the 
vertical mean flow. The results are; 
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(13) 



(14) 



In (13) and (14), p', which represents the departure of the pressure from the 
vertical average, is, using (4), expressed in terms of p as: 



0 1 f 

p'=\pgds- — ^ fpgds 



dz. 



( 15 ) 
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The method of solution consists of predicting V^y/,y/,u' ,v\T, and S from (11), 

(12), (13), (14), (6) and (7), respectively. The total current is then obtained by adding the 
vertical shear part to the vertical average part, after the latter is obtained from i}/ using (9) 

and (10). The diagnostics p, w, and p' are then obtained explicitly from the equation of 
state (5), continuity equation (8), and hydrostatic relation (15), respectively. 

C FORCmG CONDITIONS AND EXPERIMENTAL DESIGN 

Previous experiments by Batteen et al. (1992) investigated the role of steady 
equatorward winds with anticyclonic wind stress curl in the generation of features off the 
IP. In a more recent study by Batteen (1997) the role of seasonal wind forcing and 
irregular coastline geometry has been investigated in the generation of features such as 
eddies and filaments in the California Current system. Following Batteen (1997), in this 
study, seasonal wind forcing and irregular coastline geometry will be used to investigate 
the generation of similar features in the CCS. 

To explore the effects of seasonal wind forcing, the model is forced from rest with 
climatological wind fields that originate from a 2.5° by 2.5° grid of the European Centre 
for Medium-Range Weather Forecasts (ECMWF) near-surface wind analyses (Trenberth 
et al, 1990). The monthly mean stresses based on twice daily wind analyses from 1980- 
1989 have then been interpolated spatially to the 9 by 11 km model resolution and 
temporally to daily wind values. 
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An example of the wind fields used can be seen in Figure 2, which shows the 
annual migration of the Azores Subtropical High from the south in the winter (Figure 2d), 
to its most northern extent in the summer (Figure 2b). The atmospheric pressure pattern 
for November (Figure 2d) depicts a low to the north and the Azores High to the south, 
which results in a wind divergence at ~40° N. This pattern of weakly poleward winds 
north of 40° N and equatorward winds to the south continues through December. During 
January and February (Figure 2a) the divergence in the wind field migrates poleward. By 
March an equatorward component in the wind field is observed along the entire domain. 
The strongest equatorward winds are discernible from July (Figure 2b) through August. 
By September (Figure 2c) the winds start to weaken throughout the domain, and 
divergence in the wind field is observed in the north in October. In November the wind 
divergence returns to ~40° N. 

The design of the model experiments is as follows. To focus on the role of 
seasonal winds and irregular coastline geometry, experiment 1 looks at the model response 
to seasonal wind forcing along an idealized "straightened" coastline (Figure 3a), while 
experiment 2 investigates the model response to the same wind forcing as in experiment 1, 
but with an irregular, more realistic coastline (Figure 3b). In both experiments the model 
integrations start from a state of rest and, once a day, the model is updated with 
interpolated ECMWF winds. To allow the model to reach quasi-equilibrium, each model 
integration is allowed to run for 3 years. 
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m. RESULTS FROM MODEL SIMULATIONS 



A. EXPERIMENT 1 - EFFECTS OF WIND FORCING AND IDEALIZED 

COASTLINE GEOMETRY 

1. Spin-up Phase 

In response to the predominantly equatorward winds over the model domain, a 
coastal equatorward surface current develops along the entire coast by day 45. The 
current develops initially in the southern end of the domain in response to the stronger 
equatorward winds associated with the southerly position of the Azores High. By day 180 
(Figure 4a) there is generally equatorward flow throughout the domain. In addition, a 
poleward undercurrent has also developed, initially, in the southern end of the domain and 
then extending along the entire coast. A cross-section of meridional velocity (Figure 5a) 
shows the typical vertical structure of the equatorward coastal jet and the poleward 
undercurrent. The coastal jet core remains within ~50 km of the coast, extending offshore 
to -200 km. The core velocity ranges from -30-35 cm/s, and extends from -100 m depth 
near the coast to -500-600 m depth offshore. The weaker undercurrent is also seen with a 
core velocity of —5 cm/s lying beneath the coastal jet within -40 km of the coast. 

As the wind becomes equatorward over the entire domain and increases in 
intensity, increased stress is exerted on the ocean surface creating Ekman transport 
offshore. This results in upwelling of cooler water along the coast. In the spring, 
upwelling is confined to the south where the stronger winds exist, while by summer 
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upwelling is evident along the entire coast with fairly vigorous upwelling occurring just 
south of Cabo de Sao Vincent and south of Cape Beddouzza/Cape Ghir (e.g., Figure 4a). 

As the coastal jet and undercurrent become fully established (for example, see 
Figure 5b, which shows the structure of the currents just prior to meander formation), the 
currents become unstable and form meanders. This occurs off northwest Africa first 
where there has been more time for the strong equatorward winds to act on the ocean 
surface. By day 255 (Figure 4b) a meander develops just upstream of Cape Beddouzza 
and by day 270 (Figure 4c) forms an eddy inshore of the coastal jet in the cold upwelled 
water. No meanders form in the northern portion of the model until approximately day 
330 (Figure 4d). Since this meander develops in the late fall when the equatorward winds 
are diminishing (thus weakening the equatorward coastal current), this meander is not as 
extensive as that to the south off northwest Africa. 

In the fall, the currents begin to weaken with the slacking winds. The upwelling 
also ends, first in the poleward end of the domain and then progressively equatorward 
until by around day 345 all upwelling has ceased. By approximately day 336 (Figure 4e) a 
surface poleward current appears in the extreme northern portion of the model domain in 
response to the weak poleward winds present off the northwestern coast of the IP in 
November and December. 

In the spring of the second year, as the Azores High ridges back to the north, 
upwelling of cool water redevelops along the coast as the equatorward flow reestablishes 
itself and intensifies into the summer. In the second year, the coastal currents become 
unstable sooner, developing meanders earlier. To the south, the meanders are larger. 
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while to the north along the IP, the meanders become large enough to create eddies 
inshore of the coastal jet. Filaments also form as the coastal jet meanders through the 
anticyclonic eddies. The filaments then flow offshore and advect the cold upwelled water 
near the coast offshore. By April a large-scale anticyclonic circulation sets up in the Gulf 
of Cadiz and a fairly vigorous filament develops west of Cape Beddouzza. In May (Figure 
4f) a cold core cyclonic eddy forms on the offshore extension of this filament. Through 
the remainder of the year the equatorward jet, embedded with eddies, continues to 
meander downstream. 

After allowing the model to run long enough, it can be seen that the system 
reaches a quasi-equilibrium state and that many of the features continue to be generated 
and sustained. By the beginning of the third year, the mean kinetic energy has leveled off 
reaching a quasi-equilibrium state. This is quasi-equilibrium state is expected due to the 
seasonal variation of the wind forcing which is the source of the mean kinetic energy. 

2. Quasi-equilibrium Phase 

By year three, the model has reached a quasi-equilibrium state. To see the 
seasonal structures of the features of the CCS, the model output fields are time-averaged 
every month and also for the upwelling season. 

In the spring (Figure 6a), a relatively strong equatorward current of ~30-35 cm/s 
exists ~200-250 km off the coast of Iberia lying near the 17.5° C isotherm, with a stronger 
coastal jet appearing closer to shore in regions where the current bends anticyclonically 
near the coast (e.g. between 37° N and 38° N in Figure 6a). A large anticyclonic eddy 



17 



that set up in the previous year in the Gulf of Cadiz remains. Off the northwest coast of 
Africa a series of large amplitude meanders exist with several cyclonic and anticyclonic 
eddies on the order of 150 km. Upwelling begins in the spring off Cape Beddouzza, 
spreading north through late spring and into summer. As the upwelling intensifies, the 
coastal jet advects the cold water offshore, forming a distinctive filament stretching 
westward to ~250 km. 

By summer (Figure 6b), upwelling has spread north along the entire coast and in 
August reaches a maximum along the IP. A strong coastal equatorward flow also 
reappears along the entire IP in conjunction with the continued equatorward wind stress. 
By early summer the coastal jet becomes unstable, forming a meander at -38.5° N. By 
July, two distinct meanders exist, one still near -38.5° N and a second off Cabo de Sao 
Vincent, with cyclonic cold core eddies inshore of the currents at both locations. Along 
northwest Africa maximum upwelling occurs a month earlier than off the IP. Large 
meanders exist in response to the stronger currents. As the strengthened currents flow 
offshore in the regions of the developed meanders, the current advects cold upwelled 
water offshore with the result that several extensive filaments extend offshore for 
hundreds of kilometers. 

Persistent patterns for the entire year 3 upwelling season were investigated by 
averaging the temperature and current fields over that season. The average maps (Figure 
7a and 7b) reveal that an offshore equatorward meandering current, an equatorward 
coastal jet, a poleward undercurrent, eddies, and filaments, appear as regular features of 
the CCS during the upwelling season. Both the offshore current and coastal jet meander 
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cyclonically and anticyclonically. The offshore current has alongshore wavelengths of 
several hundred kilometers, cross-shore excursions of -100-200 km (with a large 
anticyclonic meander occurring to the west of the Gulf of Cadiz), speeds of 15-25 cm/s; 
and a vertical excursion to -400-600 m depth. The coastal jet flows alongshore of the 
Iberian coast until at —39° N it develops a meander with a wavelength of —150 km. Off 
Cabo de Sao Vincent, the coastal jet veers offshore again to form a second meander. The 
coastal jet then continues in a nearshore anticyclonic flow through the Gulf of Cadiz. At 
-34° N a third meander exists along the coast of northwest Africa. The coastal jet has 
speeds of -20-35 cm/s and extends to -400 m depth. The filaments extend -300 km 
offshore with widths of -30-100 km. The eddies are on the order of 100-150 km in 
diameter and tend to extend fairly deep, below 600 m depth (e. g., see Figure 7b). Below 
the coastal jet lies the poleward undercurrent with speeds of -5 cm/s (Figure 7b). 

Figures 8a and 8b are horizontal maps of the upper layer mean kinetic energy 
(MKE) and eddy kinetic energy (EKE) averaged over the upwelling season. Holland et al. 
(1983) showed that maps of MKE and EKE can be used to locate sources of the mean and 
eddy energy. Comparisons of Figures 7a and 8a show that high values of MKE follow the 
coastal jet, with lower values further offshore (at -13° W) near the core of the offshore 
equatorward flow. The highest values of MKE appear near areas of tightened temperature 
gradients (i.e. near areas of greatest upwelling) and high current velocities, with the 
maximum occurring off the Cape Beddouzza region. In contrast, a comparison of Figure 
7a and 8b, shows that the highest values of EKE are found south of -37° N where the 
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largest cyclonic meanders occur. Other areas of EKE are found near areas of weak 
cyclonic meanders off the IP and in the large-scale cyclonic flow west of the IP. The 
maximum EKE appears west of the Cape Beddouzza area in conjunction with the strong 
coastal jet turning equatorward as it merges with the equatorward flow further offshore. 
A comparison of Figures 8a and 8b shows that on average the MKE is larger than the 
EKE during the upwelling season; however, the EKE covers a larger portion of the 
domain. This is consistent with other observations of the model data, which show that 
eddies are generated from baroclinic instability due to vertical shear between the surface 
equatorward current and the poleward undercurrent, as well as from barotropic instability 
due to horizontal shear between the currents. 

By fall (Figure 6c), upwelling has weakened along the entire coast. The strong 
equatorward coastal jet weakens first, moving offshore along with the temperature front. 
By November (Figure 6c), weak poleward and equatorward flows exist out to ~350 km 
where the stronger equatorward core is located. Relatively strong equatorward flow still 
exists near the IP to the south, where a large anticyclonic meander forces the current near 
the coast. A meander off Cabo de Sao Vincent with a cold core cyclonic eddy inshore of 
the current also remains through the fall as does the large scale anticyclonic eddy in the 
Gulf of Cadiz. Off northwest Africa major meanders are maintained, as are several large 
filaments. 

By winter (not shown), only a weak equatorward flow exists over the entire 
domain. Embedded in this flow are numerous jets, eddies, and filaments. 
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On the basis of these results it is seen that a seasonal cycle occurs for the coastal 
currents, upwelling, and filaments. The eddies can be seasonal as well as quasi-permanent, 
while the meanders can be quasi-permanent. These quasi-permanent features appear to 
play a significant role in the modification of the seasonal features which leads to large 
temporal and spatial variability in the CCS. 

B. EXPERIMENT 2 - EFFECTS OF WIND FORCING AND IRREGULAR 

COASTLINE GEOMETRY 

1. Spin-up Phase 

In experiment one, an idealized, ’’straightened” coastline was used. In experiment 
2, a more realistic coastline with extensive coastline irregularities is used to investigate the 
effects of capes and bays on the formation and sustainment of the features found in the 
CCS. The same temporally and spatially varying winds as those used in experiment 1 are 
used to force the model in experiment 2. All other parameters are the same as in 
experiment 1. 

The current and temperature (e.g.. Figures 9 and 10) patterns in experiment 2 
show that many of the same features that exist in experiment 1 are also present in 
experiment 2. Both an equatorward surface current and poleward undercurrent develop. 
The equatorward flow begins in the south and then progresses north as the equatorward 
wind forcing migrates north with the Azores High. In addition, as the core of the 
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undercurrent intensifies in conjunction with the progression of the upwelling season, it 
shallows and displaces the core of the equatorward surface current further offshore 
(Figures 10a and 10b) as in experiment 1 (e.g.. Figures 5a and 5b). Strong vertical and 
horizontal shears in the upper layers develop between the opposing currents, causing the 
flow to become both barotropically and baroclinically unstable, which leads to the 
formation of meanders, eddies, and filaments. 

As in the spin-up phase in experiment 1, an equatorward surface current develops 
rather uniformly along the entire coast. While the upwelling in experiment 1 was more 
uniform, only having increases near the two prominent capes (e.g., Cabo de Sao Vincent 
and Cape Beddouzza) depicted in the idealized coastline, in experiment 2 the 
corresponding upwelling is quite patchy, appearing most intense near even the less 
prominent capes (e.g., off Cabo de Roca in Figure 9a and 9b). Since these capes are areas 
where the alongshore component of the wind stress is at a local maximum, the coastal 
current velocity is also at a maximum, which leads to enhanced upwelling (Batteen, 1997). 

In experiment 1 the majority of features had no geographical preference other than 
the features that appeared tied to the major geographical features depicted in the idealized 
coastline (e.g., the cyclonic eddy off Cabo de Sao Vincent, the anticyclonic flow through 
the Gulf of Cadiz, and the large filament off Cape Beddouzza). In contrast, experiment 2 
frequently features additional eddies and filaments “anchored” off of less prominent 
coastal features (e.g., off Cabo da Roca, Cabo de Sines, Cabo de Santa Maria, and the 
Cape Beddouzza/Cape Ghir area shown in Figure 9d) that occur along the realistic 
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coastline. In addition, in experiment 2 upwelling occurs significantly sooner (~ two 
weeks) and with greater intensity off the DP, where the largest differences occur between 
the idealized and realistic coasts (e.g., compare Figures 4 and 9). Experiment 2 also 
develops more meanders in the vicinity of the additional capes of the realistic coastline 
which also form earlier then those in experiment 1. These meanders intensify and develop 
into predominantly cyclonic cold core eddies, which in time coalesce with other eddies to 
form larger cyclonic eddies on the order of -100-300 km. Similar to that in experiment 1, 
by the end of the upwelling season (Figure 9c), the equatorward surface flow has taken 
the form of a meandering jet embedded with several eddies. - 



2. Quasi-equilibrium Phase 

Similar to experiment 1, by year 3 the model reaches quasi-equilibrium and shows 
that most features continue to be generated and maintained. As in experiment 1, to see 
the seasonal structure of the CCS, the model output is time averaged every month as well 
as during the upwelling season. 

The model results show that in the spring (Figure 11a) there is a mean 
equatorward flow embedded with numerous eddies. The core of the jet appears to follow 
the temperature gradient of 16.5-17.5° C. This jet flows close to the coast off northern 
Iberia, but leaves the coast as it impinges on Cabo de Roca. It then continues to meander 
equatorward following the temperature front south. By April (Figure 11a), a coastal jet 
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appears off the coast of Northwest Africa and, as in experiment 1, upwelling is also 
reestablished by this time. By early summer (Figure lib) upwelling begins to appear in 
pockets along the Iberian coast and a cold core cyclonic eddy develops southeast of Cabo 
da Sao Vincent. A relatively robust coastal jet (of ~20-40 cm/s) extending to depths of ~ 
600 m (Figure Ilf) is also now found along the entire coastline. It appears to flow along 
the newly established tightened temperature gradient of 14.5-15.5° C as additional cold 
water is upwelled near the coast. The maximum velocity variations are normally found 
near the capes, with increases poleward and decreases equatorward of the capes. Note the 
cooler water inshore and warmer water offshore of the coastal jet and the numerous 
eddies extending beyond 800 m depth (Figure 1 If). 

Maximum upwelling occurs by late summer (Figure 1 Ic). The strong equatorward 
coastal jets which have developed get deflected offshore by the large capes (e.g., Cabo da 
Roca) and advect the cold upwelled water offshore producing filaments. In addition to 
occurring off Cabo da Roca, filament generation is also observed off Cabo de Sao 
Vincent, Cabo de Santa Maria, and Cape Beddouzza/Cape Ghir. 

To better observe the features of the upwelling season the temperature and current 
fields are averaged over the entire upwelling season. This shows that an equatorward 
meandering current offshore, an equatorward coastal jet, a poleward undercurrent, eddies, 
and filaments appear as regular features of the CCS during the upwelling season (Figures 
12a and 12b). The surface equatorward current flows along the coast of northern IP until 
~40° N where it impinges on Cabo de Roca splitting into two distinct flows, a coastal jet 



24 



and an offshore current. Both the offshore current and coastal jet meander cyclonically 
and anticyclonically. The offshore current has alongshore wavelengths of several hundred 
kilometers, cross-shore excursions of -100-200 km, speeds of 15-35 cm/s, a vertical 
excursion to -600-800 m depth, and seems to flow approximately along the 17® C 
isotherm. The coastal jet continues to flow alongshore of the Iberian coast continuing 
anticyclonically through the Gulf of Cadiz and down the coast of Africa. The coastal jet 
has speeds of —25-40 cm/s with increases near areas of maximum upwelling and extends to 
-1000 m depth. The filaments extend -200 km offshore with widths of -30-100 km. The 
eddies are on the order of 100-150 km in diameter and tend to extend below 1000 m 
depth. Below the coastal jet lies the poleward undercurrent which is present along the 
entire coast (e.g.. Figure 12c) throughout the upwelling season. The undercurrent 
meanders within -100 km of the coast at speeds of -10-15 cm/s. Offshore of the 
undercurrent are eddies of 0(100 km) diameter. 

Figures 13a and 13b are horizontal maps of the upper layer MKE and EKE 
averaged over the upwelling season. As in experiment 1, MKE is found principally near 
the core of the surface currents, while EKE is preferentially found in areas of cyclonic 
meanders and eddies. In addition, MKE has larger values than EKE; however, the EKE 
extends over a larger region. Comparisons of Figures 8a and 13a show that the MKE is 
spread out over a larger region than in experiment 1, while the maximum MKE has shifted 
from northwest Africa to the IP. Comparing Figures 8b and 13b shows that the large 
region of EKE west of the Gulf of Cadiz in experiment 1 has shifted north in experiment 2 
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to west of the IP. In addition, a large region of EKE is found in experiment 2 along the 
northern portion of the IP to Cabo da Roca as well as off the southern coast of the IP. 
Large areas of EKE south of the Cape Beddouzza area are found in both experiments. 
The regions with the largest values of EKE can be interpreted as areas where eddies are 
likely to be generated during the upwelling season. 

In the fall (Figure lid), the coastal equatorward flow moves offshore of Iberia as 
the temperature front weakens and is displaced offshore. The eddies that formed inshore 
of the current earlier also begin to propagate westward. Further offshore, a stronger 
equatorward flow is present (-30-40 cm/s) embedded with numerous cyclonic and 
anticyclonic eddies. To the south off Cape Beddouzza and Cape Ghir a strong filament 
forms with a cold core eddy at its western extent and a warm core anticyclonic eddy to the 
northeast. 

In the winter (Figure lie) the coastal jet is displaced even farther offshore by the 
continuing westward propagation of the nearshore eddies. By November a poleward flow 
(-20 cm/s) is established in the wake of the westward flow of the surface current and 
eddies (Figure 1 Ig). This flow lasts through December and is likely due to forcing by the 
weak poleward winds of November (Figure 2d) and December. 

Experiment 2 shows that irregularities in the coastal geometry are key elements for 
"anchoring" upwelling and filaments to geographical features. It also shows that coastline 
irregularities aid in the development and growth of meanders and eddies. 
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c. 



COMPARISONS OF MODEL RESULTS WITH OBSERVATIONS 



It is useful to do a comparison of model results with available field and satellite 
observations to ensure they are consistent with actual features of the CCS. Although 
direct comparisons are not valid since the studies are not hindcasts, but idealized process- 
oriented studies, it can be shown that the phenomenonological model behavior is 
qualitatively similar to the observed data. General features of the Canary Current are 
apparent from satellite images as well as measured current data. The discussion in the 
ensuing paragraphs details the general features, and then discusses the specific magnitudes 
and characteristics along the western coasts of the IP and Northwest Africa. Overall, the 
results of the model simulations of experiment 2 highlight the major characteristics of the 
CCS with impressive similarities to field observations. 

1. Comparison of Ocean Currents 

The major characteristics of the modeled Canary Current off western Iberia and 
northwestern Africa are consistent with the few field observations and studies available. 
Satellite imagery also reveals a system of meanders, eddies, and along-shelf jets, on a 
range of length scales that correspond closely with model output. 

a. Equatonvard Flow 

In experiment 2, the numerical model, through the use of climatological 
wind data and a realistic coastline, develops an ocean current that contains many features 
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observed in the CCS. The existence of a mean broad equatorward flow throughout the 
entire model domain matches the general description of classical EBCs by Wooster and 
Reid (1963). In addition, the sustained magnitudes of the model’s equatorward flow 
compares favorably with available data. For example, Figures 14a and 14b from 
DMAH/TC (1988) show a mean flow of 0.5 to 0.7 kts (~25-35 cm/s) with occasional jets 
reaching 2 to 2.5 kts (~100-125 cm/s), which are consistent with velocities found in the 
model simulations. 

b. Poleward Flow 

The model establishes a poleward undercurrent with a core at ~300-400 m 
depth but extending to over 1200 m depth, with magnitudes consistent with observations 
from current measurements noted by Haynes and Barton (1990) of ~20 cm/s with a core 
depth ranging between 100 and 600 m. Ambar (1982) discusses the existence of three 
separate cores (~400 m, 600-900 m, and 1100-1200 m depth) of poleward MO (with 
velocities of ~18-29 cm/s) existing along the southern coast of the IP in the Gulf of Cadiz 
and extending to the north along the entire western coast of the IP. The model depicts 
this deep poleward undercurrent (-10-15 cm/s) as well (Figures 5, 7b, 12b, and 12c). 

The model also depicts a narrow (-50 km), weak (-15-20 cm/s) surface 
poleward flow, extending to -400 m depth. The flow is found close to shore in the 
northern portion of the model domain during December (Figures lie and llg). This 
corresponds well to the shift in the climatological winds from northerly to southwesterly 
allowing the equatorward flow to diminish. Recent studies by Frouin et ah, (1990) and 
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Haynes and Barton (1990) using high resolution infra-red images from NOAA7 and 
NOAA9 and numerous in-situ measurements reveal the existence of a surface poleward 
flow off the northern coast of Portugal and northwestern Spain. These studies show a 
current -25-40 km wide, with speeds ~25 cm/s, extending to ~200 m depth. 

2 . Comparison of Eddies 

Observations by Fiuza (1984) show pairs of mesoscale eddies on the order of 100 
km exist off the west coast of Iberia (Figure 15). This study, as well as new high 
resolution satellite imagery confirms the models generation of numerous eddies in the 
CCS. 

In experiment 2, the onset of instabilities occurs at preferred locations. These 
locations are where the coastline significantly changes, i.e., near Porto, Cabo da Roca, 
Cabo de Sines, Cabo de Sao Vincent, the Gulf of Cadiz, and Cape Beddouzza/Cape Ghir 
(Figure 1). Figure 15 shows a satellite image that agrees with some of the west coast 
preferred locations for the formation of meanders and eddy positions, while Figure 16 
shows that the model also forms eddies in similar locations. 

3. Comparison of Upwelling 

Because the currents are dependent on the winds and the winds have been shown 
to oscillate seasonally, upwelling off the coast of Portugal (Fiuza et ai, 1982) as well as 
along the coast of the rest of the model domain also varies with the seasons. Upwelling is 
observed off Portugal from -43° N to Cabo de Sao Vincent in the south (Figure 17) and 
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along the southern coast of Portugal to Cabo de Santa Maria, ~8° W (Figure 1 8), which 
the model depicts quite well (e.g.. Figure 1 lb). Observations of upwelling at its maximum 
during the summer and early autumn are shown in Figures 19 and 20, while Figure 21 
shows the model depicting increased upwelling during this period as well. Weak to non- 
existent upwelling is observed to occur along the northern coast of Morocco, then 
increases again south of 33° N, near Cape Beddouzza (Figure 22)(Fiuza 1982). As the 
Azores High migrates north, the upwelling first appears off Cape Ghir/Cape Beddouzza in 
early May as seen in Figures 23a and 23b of Van Camp et al. (1991), then gradually 
moves north up the Moroccan coast. This is simulated by the model and depicted in 
Figures 11b, 21, and 24. 

4. Comparison of Filaments 

Several studies have been done that show the existence of cold filaments extending 
off the coasts depicted in our model: Fiuza (1983), Barton (1986), and Fiuza and Sousa 
(1989), off the northwest coast of Iberia; Folkard et al, (1997) off the southern coasts of 
Portugal and Spain in the Gulf of Cadiz; and Van Camp et al, (1991) and Hagen et al, 
(1996) near Cape Ghir off northwest Africa. 

a. West Coast of Iberia 

Fiuza (1983) was the first study indicating the existence of filaments off the 
west coast of Iberia. Additional studies by Fiuza and Sousa (1989) confirmed the 
existence of filaments using CTD data and satellite imagery in which several filaments 
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appear including a large filament off Cabo da Roca (Figure 25). In addition, recent sea 
surface temperature (SST) observations show several good images of filaments occurring 
along the entire western coast of Iberia (Figure 19) with a particularly long filament near 
40.5° N. The model depicts filaments in many of the same locations as well (Figure 21). 
A rather strong filament exists near -40.5° N as well as one further south at -39.5° N. A 
small filament is also evident at -37.5° N. Figure 11a shows a filament off Cabo da Roca 
with a striking resemblance to the filament in Figure 25. 

b. South Coast of Iheria/Gulf of Cadiz 

Folkard et ai, (1997) investigated the waters of southern Iberia in the Gulf 
of Cadiz. In that study he observed filaments extending off the coast during March and 
October (Figures 26 and 27). The model depicts these filaments as well: a coastal jet 
deflected offshore near -8° W, Cabo de Santa Maria, advects cold upwelled water 
offshore (Figures 28 and 29). 

c. Cape Ghir/Northwest Africa 

Consistent with satellite observations (Figure 23) of Van Camp et al, 
(1991) and satellite observations (Figure 30) coupled with CTD surveys (Figure 31) from 
Hagen et ai, (1996), a filament of cold upwelled water extends off Cape Ghir near 3 1° N 
for about 200 km. The model depicts this filament as an equatorward meandering coastal 
jet which transports cold upwelled water offshore. Observations show that the coastal jet 
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and subsequent frontal zone caused by the advected cold water results in the generation of 
a cyclonic eddy-like feature on the cold water side, while an anticyclonic eddy forms on 
the warm water side (Figure 32). Most of the general features of the filament are depicted 
well in the model; however, the placement and orientation of the filament are slightly off. 
The model places the filament extending offshore at about 32° N with a 
northwest/southeast orientation. This discrepancy may be due to the lack of bottom 
topography in the model since topographic influence of the Cape Ghir Plateau on the 
actual formation of the filament has been demonstrated by Hagen et al. (1996). 

The variety of length scales and life time scales, in conjunction with the 
differences in field observations themselves, emphasize the complex flow region of the 
CCS. Except for the slight disagreement with placement and orientation of the cold 
filament off Cape Ghir, the numerical model results of experiment 2 show striking 
similarities to the field observations. 
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IV. SUMMARY 



This study was designed to investigate the role of wind forcing and coastal 
geometry in the formation of currents, eddies, and filaments in the Canary Current system. 
It was found that wind forcing can set up coastal currents, which can subsequently become 
unstable and lead to the formation of many of the observed features in the CCS. To 
demonstrate that wind forcing can create both a surface current and an opposing 
undercurrent, a high resolution, multi-level, PE model was forced from rest with spatially 
and temporally varying winds. The major components of the currents developed quite 
rapidly, generating both horizontal (barotropic instability) and vertical (baroclinic 
instability) shears which caused the currents to become unstable, resulting in the formation 
of meanders. As the meanders grew, eddies developed along with areas of increased 
upwelling. In time, filaments also formed as the meanders intensified and advected cold 
upwelled water offshore. 

A seasonal cycle was also identified, as was expected due to the temporal nature of 
the seasonal winds. Obviously there is a distinct upwelling season tied to the migration of 
the Azores High. In addition to this, however, is the generation of a seasonal surface 
equatorward coastal jet overlying a poleward undercurrent that occurs during the 
upwelling season. The eddy formation along the coast also appears to be seasonal, 
forming during the height of the upwelling season, when the coastal jet is at its maximum 
strength. In addition, filaments appear to be seasonal, occurring during and just after the 
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upwelling season; however, some filaments do occur during other seasons, especially off 
Iberia where cold water exists nearshore throughout the year. 

The meanders and some eddies appear as quasi-permanent features. A relatively 
large anticyclonic eddy develops in the Gulf of Cadiz and remains there year-round. 
Additional quasi-permanent meanders and eddies play a significant role in the modification 
of the coastal current, upwelling, and filament formation. 

Irregularities in the coastal geometry also appear to be significant factors in the 
formation of certain features as well as for enhancing the growth of meanders, eddies and 
filaments. Specifically, capes are significant in the location and formation of filaments, 
providing an area of enhanced upwelling as well as acting a barrier to force the coastal jet 
offshore. 

A comparison of model results with available observations show that while both 
experiments successfully simulate most of the features of the CCS, experiment 2 shows 
better agreement between current velocities and feature placements. Although experiment 
2 simulated most features well, some discrepancies existed. For example, the positioning 
of the Cape Ghir Filament placement in the model run is approximately half a degree north 
of its actual position. This may be due to the lack of bottom topography in the model, 
since Hagen et al. (1996) suggests that the Cape Ghir Plateau influences the formation of 
the filament. Another example of placement error due to lack of bottom topography is the 
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placement of the poleward surface current off northern Iberia. The model places the 
current approximately one degree closer to the coast than available observations show. 

Besides the lack of bottom topography the model also does not include the effects 
of salinity forcing or an opening of the boundary to the Mediterranean Sea. Several 
variable salinity inputs exist in the area, including the Mediterranean Sea influence, as well 
as several large river inflows off northern Iberia and numerous smaller rivers along 
Morocco. Opening the boundary to the Mediterranean Sea would allow the existence of 
MO, which would then result in a density-driven component to the undercurrent (in 
addition to the model-simulated wind driven component). In addition, incorporation of 
bottom topography, salinity forcing, and MO could result in the proper simulation of other 
cores of the undercurrent, which have been observed in the region (e. g., Ambar, 1982). 

While this study focused on the role of seasonal wind forcing and used the 
relatively coarse spatial resolution of ECMWF winds, it would be of interest to investigate 
the role of higher resolution spatial and temporal wind forcing. Fine spatial resolution 
down to ~1 km combined with fine temporal resolution of ~1 day would allow the oceanic 
response to wind events, relaxations, and reversals to be studied. Observations have 
shown that, depending on the prevailing wind direction, an inflow or outflow through the 
Strait of Gibraltar can exist, which could affect the flow in the Gulf of Cadiz (e.g., Folkard 
etal, 1997). 
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Of concern to studies such as this one is the lack of readily available data along 
both the Moroccan coast and the IP. Few observations exist of the equatorward flow, 
including velocities and structure. Some data exists on the relatively deep poleward MO 
and the upwelling along Iberia, and more recently several studies have discussed the 
poleward surface current located off northern Iberia. In addition, several studies discuss 
briefly the surface temperature, salinity, and density patterns off northwest Africa and of 
the Cape Ghir Filament, but little discussion is available about the basics. This area is a 
prime candidate for a refreshing update of the generalization given the Canary Current by 
Wooster and Reed (1963). 

Nevertheless, the results from these model runs show surprising similarities to the 
large-scale observational data available in the area studied. This supports the hypothesis 
that wind forcing coupled with irregular coastal geometry are key issues in the generation, 
evolution, and sustainment of the currents, and other unique feature including meanders, 
eddies, and filament structures of the CSS. 
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Figure 1. Domain of the model for the Canary Current System (CCS) off the western 
coast of the Iberian peninsula and northwest Africa. The domain is bounded by 30° N to 
42.5° N, 5° W to 17.5° W. The figure also shows the major geographical features. 
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Figure 2. Climatological (1980-1989) ECMWF winds in m/s used in experiments 1 and 2: 
(a) February, (b) July, (c) September, (d) November. 
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Figure 3. Coastlines used in experiments 1 and 2: (a) The idealized or "straightened" 
coastline uses a straight coast for the Iberian Peninsula with a 90 corner depicting Cabo de 
Sao Vincent. The Gulf of Cadiz is then depicted with a straight coast for southern Iberia 
as well as a portion near the Strait of Gibraltar, and then a sloping straight coast south to 
an angle depicting the Cape Beddouzza area, (b) The irregular or "realistic" coastline 
which features the general characteristics and irregularities of the coast in the model 
domain as described in Figure 1 . 
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Figure 4. Experiment 1 : Temperature contours and velocity vectors at 13 m depth at 
days (a) 180, (b) 255, (c) 270, (d) 330, (e) 336, and (0 501 (day 136 of year 2). The 
contour interval is 1°C. To avoid clutter, the velocity vectors are plotted every third grid 
point in the cross-shore direction and every fourth grid point in the along-shore direction. 
Maximum current velocity is 100 cm/s. 
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Figure 5. Experiment 1; Cross-shore section at 38,8° N of the meridional component of 
velocity (v) in the coastal region at days (a) 255 and (b) 309. Solid lines indicate poleward 
flow, while dashed lines indicate equatorward flow. The contour interval is 2.5 cm/s for 
poleward flow and 5 cm/s for equatorward flow. 
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Figure 6. Experiment 1: Temperature contours and velocity vectors at 13 m depth in the 
third year of the model simulation time-averaged over the months of (a) April, (b) August, 
and (c) November. The contour interval is 1° C. Maximum current velocity is 100 cm/s. 
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Figure 7. Experiment 1; Time-averaged plots for the upwelling season (April-September) 
of year 3 of (a) temperature contours and velocity vectors at 13 m depth, and (b) cross- 
shore section of velocity (v) at 37.7° N. The contour interval is 1° C in (a) and 2.5 cm/s 
for poleward flow and 5 cm/s for equatorward flow in (b). Maximum current velocity is 
100 cm/s. 
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Figure 8. Experiment 1; (a) Mean and (b) eddy kinetic energy at 13 m depth. Contour 
interval is 100 cm^/s^. 
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Figure 9. Experiment 2: Temperature contours and velocity vectors at 13 m depth at 
days (a) 171, (b) 240, (c) 300, and (d) 360. The contour interval is 1° C. Maximum 
current velocity is 100 cm/s. 
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Figure 10. Experiment 2: Cross-shore section at 39.5° N of the meridional component of 
velocity (v) in the coastal region at days (a) 255 and (b) 309. Solid lines indicate poleward 
flow, while dashed lines indicate equatorward flow. The contour interval is 2.5 cm/s for 
poleward flow and 5 cm/s for equatorward flow. 
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Figure 11. Experiment 2: Temperature contours and velocity vectors at 13 m depth in the 
third year of the model simulation time-averaged over the months of (a) April, (b) June, 
(c) August, (d) October, (e) December, and cross-shore section of velocity (v) at (f) 37.8° 
N time-averaged over June and (g) 41.1° N time-averaged over December. The contour 
interval is 1° C in (a) - (e), 2.5 cm/s for poleward flow and 5 cm/s for equatorward flow in 
(f) and 5 cm/s for (g). Maximum current velocity is 100 cm/s. 
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Figure 12. Experiment 2: Time-averaged plots for the upwelling season (April- 
September) of (a) temperature contours and velocity vectors at 13 m depth, (b) cross- 
shore section of velocity (v) at 32.1° N, and (c) temperature contours and velocity vectors 
at 870 m depth. The contour interval is 1° C in (a) and 2.5 cm/s for poleward flow and 5 
cm/s for equatorward flow in (b). Maximum current velocity is 100 cm/s in (a) and 50 
cm/s in (c). 
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Figure 13. Experiment 2; (a) Mean and (b) eddy kinetic energy at 13 m depth. Contour 
interval is 100 cmW. 



75 



DEPTH 



13m 



DATA SET: meonupw_yr3 




41 0"N 



39 0"N 



37 

bJ 

Q 

t 

h- 

< 



35 C*N 



33 0*N 



3i 0*N 



17 0“W 



15 0*W 



13 0“W 



9 C^W 



7 0”W 



5 0^W 



CONTOUR IVASK 



LONGITUDE 

EDDY KINETIC ENERGY ((cm/sec)**2) 



b 



76 




Figure 14. Mean surface current speeds in knots and prevailing direction for (a) winter 
(January, February, March), for (b) summer (July, August, September). "After 
DMAH/TC, 1988." 



77 








' ' t / ^ ftc >0 fT 

■ /:■ - : 



/• 

& 



SfE PlO" t T^ 



V^\ 

i 





LEGEND 


ARROW POSITION INDICATES THE DIRECTION TO 

WARD which the prevaiung surface currents 
FLOW (SET). WHILE its COLOR INDICATES THE 
percentage of time THE PREVAILING SURFACE 
CURRENTS FLOW THE SPECIFIED DIRECTION AND 
SPEED (CONSTANCY OF THE CURRENT). 

FIGURES FOR SPECIFIED AREAS ARE DISPLAYED ON 
THE FOLLOWING PAGES. 


— 


<20% (VARIABLE) 




/VTO 1 V 7b 

A0% TO 55% 

56% TO 75% 

>75% 

REGIONS OF SIGNIFICANT TIDAL CUR 
RENTS 


— 


□ 


0 s 


MEAN SPEED (KNOTS) OF THE PREVAILING 
SURFACE CURRENTS, 



>>?.MxLU^iJ.k. 
\V0"E 10" 20^ 



.■■liJ..li.:.LUl I ! I : /i5 

30" 40" 



b 



78 




Latitude 



40° O’N 



Eddy Formation Locations 




Figure 15. Computer representation of AVHRR NOAA7 Sea Surface Temperatures 
(SST) with observed eddy formations from 5 August 1982. "After Fiuza 1984." 
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Figure 16. Experiment 2: Temperature contours and velocity vectors at 13 m depth at 
day 984 showing preferred eddy formation locations. The contour interval is 1° C. 
Maximum current velocity is 100 cm/s 
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Figure 17, AVHEUR. SST from 21 June 1995 showing increased cold SST southwest of 
Cabo de Sao Vincent caused by intensified upwelling and advection of upwelled water 
around the cape. "From The Remote Sensing Data Analysis Service (RSDAS), Plymouth 
Marine Laboratory, UK." 
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Figure 18. AVHRR thermal image from 7 July 1989, showing effects of westerly winds 
on the surface temperature fields in southern Iberian coastal waters during the summer. 
"From Folkard et al., 1997." 
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Figure 19. AVHRR SST from 27 July 1993 showing strong upwelling and numerous 
filaments off the west coast of the IP. "From The Remote Sensing Data Analysis Service 
(RSDAS), Plymouth Marine Laboratory, UK." 
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Figure 20. Sea Surface Temperatures (SST) showing intense upwelling north of Cabo da 
Roca from 29 July 1977. "From Fiuza 1982.” 
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Figure 21. Experiment 2: Temperature contours and velocity vectors at 13 m depth in the 
third year of the mode! simulation time-averaged over the month of July shows the 
increased upwelling during the summer. The contour interval is 1° C. Maximum current 
velocity is 100 cm/s. 
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Figure 22. Seasonal cycle of thermal anomalies along the western coasts of Portugal and 
northwestern Africa between 32° N and 44° N. The thermal anomalies are computed 
from the climatological sea surface temperatures of Crowe (1949) and from observations 
from the period 1959-1969. The contour interval is 0.5° C. "After Fiuza 1982." 
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Figure 23. SST from (a) 3 1 May 1987 and (b) 9 June 1987 showing strong upwelling 
progressing north up the western coast of Northwest Afnce and the development of the 
Cape Ghir filament. "From Van Camp eta/., 1991." 
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Figure 24. Experiment 2; Temperature contours and velocity vectors at 13 m depth in the 
third year of the model simulation time-averaged over the month of May shows the 
upwelling beginning to the south off Cape Ghir/Cape Beddouzza in the spring. The 
contour interval is 1° C, Maximum current velocity is 100 cm/s. 
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Figure 25. Satellite image of the SST off the west coast of Portugal and Spain, at 1611 
UT 5 September from the University of Dundee ground receiving station. White tones 
nearshore correspond to cold, upwelled water with a distinct filament extending off Cabo 
da Roca. "From Fiuza and Sousa 1989," 
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Figure 26. AVHRR image from 7 March 1988, showing a typical surface thermal pattern 
for spring. "From Folkard e/ a/., 1997." 
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Figure 27. AVHRR image from 8 October 1984, showing a typical surface thermal 
pattern for autumn. "From Folkard et ai, 1997." 
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Figure 28. Experiment 2: Temperature contours and velocity vectors at 13 m depth in the 
third year of the model simulation time-averaged over the month of March shows the 
filament developing southwest of Cabo de Santa Maria. The contour interval is 0.2° C. 
Maximum current velocity is 1 00 cm/s. 
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Figure 29. Experiment 2; Temperature contours and velocity vectors at 13 m depth in the 
third year of the model simulation time-averaged over the months of September and 
October shows the filament originating from the cold upwelled water southwest of Cabo 
de Sao Vincent and extending to the south of Cabo de Santa Maria. The contour interval 
is 0.5° C. Maximum current velocity is 100 cm/s. 
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Figure 30. Rough contour of the bottom topography and orography in the northeast 
Atlantic Ocean; the monthly averaged SST-map (5 September - 5 October, 1992) is 
superimposed; a mean signal of coastal upwelling occurs off Portugal and Morocco; 
Hydrographic measurements are carried out along the zonal section at 32° N (15-18 
September, 1992) and within the box (29 September 3 October 1992) covering the 
extended cold-water filament above the Cape Ghir Plateau (CGP) at about 3 1° N; note 
that the plateau is a submarine extension of the Atlas Chain onland. "From Hagen et al., 
1996." 
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Figure 31. Area of investigation (box) with the sketch for detected (full line) and Active 
(broken line) patterns of geostrophic motions within the superficial layers; arrows indicate 
the cyclonic course of the southward coastal jet and the westerward filament current (F) 
over the rough bottom topography of the CGP; isobaths in metres; shaded areas indicate 
cold water zones with an active upwelling cetre in the northeast corner; two oppositely 
rotating eddies propel upwelled water out to the open Atlantic along the filament axis at 
about 3 1° N; in addition the cyclonic eddy (CE) generates upward-directed vertical 
velocities (+) while the anticyclonic eddy (ACE) produces downward vertical 
velocities (-); data are based on CTD measurements between 29 September and 3 
October. "From Hagen e/ a/., 1996." 
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Figure 32. Experiment 2; Temperature contours and velocity vectors at 13 m depth in the 
coastal region in the third year of the model simulation time-averaged over the month of 
September shows the filament and associated eddies off Cape Ghir. The contour interval 
is 1° C. Maximum current velocity is 100 cm/s. 
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Table 1. Values of Constants Used in the Model 



Constant 


Value 


Definition 


To 


278.2°K 


Constant Reference Temperature 


^0 


34.7 


Constant Reference Salinity 


Po 


1.0276 gm cm^ 


Density of Sea Water At Tq and 6'o 


a 


2.4 X lO"* (°K)‘‘ 


Thermal Expansion Coefficient 




7.5 X 10*^ 


Saline Expansion Coefficient 


K 


10 


Number of Levels In Vertical 


Ax 


9.0 X 10^ cm 


Cross-Shore Grid Spacing 


A;. 


1.1 X 10® cm 


Alongshore Grid Spacing 


H 


4.5 X 10® cm 


Total Ocean Depth 


At 


800 s 


Time Step 


fo 


0.86 X lO'* s‘‘ 


Mean Coriolis Parameter 


g 


980 cm 


Acceleration of Gravity 




2 X 10*’cm^s‘‘ 


Biharmonic Momentum Diffusion Coefficient 




2 X 10*’ cm'' s'* 


Biharmonic Heat Diffusion Coefficient 


Km 


0.5 cm’ s'* 


Vertical Eddy Viscosity 


Kh 


0.5 cm’ s' 


Vertical Eddy Conductivity 
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